Quantum phases of a two-dimensional dipolar Fermi gas 
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We examine the superfluid and collapse instabilities ol a quasi two-dimensional gas of dipo- 
lar fermions aligned by an orientable external field. It is shown that the interplay between the 
anisotropy of the dipole-dipole interaction, the geometry of the system, and the p-wave symmetry 
of the superfluid order parameter means that the effective interaction for pairing can be made very 
large without the system collapsing. This leads to a broad region in the phase diagram where the 
system forms a stable superfluid. Analyzing the superfluid transition at finite temperatures, we 
calculate the Berezinskii-Kosterlitz-Thouless temperature as a function of the dipole angle. 
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Trapped ultracold gases are increasingly being used 
to simulate solid-state systems, where clear experimental 
signatures of theoretical predictions are often lacking 1] . 
A limitation of these gases is that the interactions are 
typically s-wave [3] whereas order parameters in solid 
state systems exhibit richer p- and d-wave symmetries. 
Recent progress in the production of trapped, cold dipo- 
lar gases fs', '3, ^ promises to change this since the dipole- 
dipole interaction is long-range and anisotropic. Impor- 
tantly, the interaction in fcrmionic hcteronuclear Fesh- 
bach molecules can be large, with electric dipole moments 
on the order of a Debye 0, H, 0] ■ This opens the door 
to experimentally reaching the superfluid phase. One 
complicating feature is that the gases tend to become 
unstable when interactions are sufficiently strong and at- 
tractive. 

In this letter, we propose confining a gas of fermionic 
dipoles of mass m in the 2D plane by a harmonic trap- 
ping potential Vz{z) = mw|z^/2. For hujz S> where 
ep = kp/2m is the Fermi energy, determined by the 2D 
density = fc|-/47r, the system is effectively 2D. The 
dipoles are aligned by an external field E (see Fig.[T]), sub- 
tending an angle Q with respect to the x-y plane. The 
interaction between two dipole moments d separated by 
r is given by V{r) = D'^{1 — 3cos^ 0rd)/r^- Here, = 
(P/Aireo for electric dipoles and drd is the angle between 
r and d. The strength of the interaction is parametrized 
by the dimensionless ratio g = 2D'^kp/{3Trep). Since 
we consider identical fermions at low temperatures T, 
additional short range interactions are suppressed. We 
show that in this configuration, the effective interaction 
for pairing has p-wave symmetry and can be made large 
without the system collapsing. This makes the 2D system 
of dipoles a promising candidate to study quantum phase 
transitions and pairing with unconventional symmetry. 

A key point is that the effective dipole-dipole inter- 
action in the x-y plane can be tuned by changing the 
angle 0. This gives rise to several interesting effects. 
For O = 7r/2, the interaction is repulsive and isotropic 
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FIG. 1: In the proposed experimental setup, aligned Fermi 
dipoles are confined in the x-y plane. The dipoles form an 
angle O with respect to the x-y plane. As the angle Q is 
reduced, the region where the interaction between dipoles is 
attractive increases. 



in the x-y plane and the dipoles are predicted to un- 
dergo a quantum phase transition to a crystalline phase 
for g ~ 27 at T = Q. As 6 is decreased, the in- 
teraction becomes anisotropic in the x-y plane and for 
< arccos(l/V3), an attractive sliver appears along the 
x axis as illustrated in Fig. H]). This gives rise to two 
competing phenomena: superfluidity and collapse [1]. 
The resulting phase diagram is shown in Fig. and 
is the main focus of this letter. 

We first examine the critical angle 0c below which the 
attractive part of the interaction overcomes the Fermi 
pressure and the gas is unstable towards collapse. The 
instability is identified with a negative value of the inverse 



compressibility k ^ = n\j-,d^£/dr 
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stability of the dipolar Fermi gas using the normal phase 
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FIG. 2: The T = phase diagram. Phases are shown in 
terms of the interaction g. For weak coupling, g < 0.725, 
the Fermi pressure stabihzes the system for all angles. For 
stronger coupling, there is a critical angle Oc below which the 
attractive interaction overcomes the Fermi pressure and the 
gas collapses. The system is superfluid and stable for Oc < 
Q < Qs- The superfluid and normal regions are separated 
by a quantum phase transition. For comparison, we also plot 
as dashed lines the phase diagram calculated with no Fermi 
surface deformation effects {a — 1). 



energy density = ^kin + ^^dir + ^cx with 



^kin 
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the kinetic energy density. We allow for the possibil- 
ity that the Fermi surface can deform |9|] due to the 
anisotropy of the interaction, by writing the T = Fermi 
distribution function as /k = 6(fc|. — fc^/a^ — fc^a^). The 
direct and exchange energies are (r = ri — 

Sdh- = ^/ ci^rirf^r2n(ri,ri)l/(r)n(r2,r2), (2) 

and 

fex = -^ y"d'rid3r27i(ri,r2)l^(r)n(r2,ri) (3) 

with n(ri,r2) = {il}^ri)ip{r2)) , -0 the dipolar field oper- 
ator, and the size of the gas. In the superfluid phase, 
there is a contribution to the energy arising from the con- 
densation energy £c oc A^/ep which is small for Aq ^ ef- 
We ignore this contribution in the following. 

In the 2D limit huj^ ^ ^f, all dipoles reside in the 
lowest harmonic oscillator level ^{z) in the z-direction. 
We then have n(ri,r2) — $(zi)$(z2)?^2d(p) with 



't-2d(p) = 
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Here, Ji{x) is the Bessel function of first order, p = 
{x,y), and f(p) = ^ joP- -f o?y^ . Evaluating the in- 
tegrals in ([2]) and ([3]), we find 
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(6) 



is a dimcnsionlcss function of the Fermi surface defor- 
mation parameter a and the dipole angle 0, and t = 

SIT? (j) + cos^ (jj/a'^. Equations (P), and ^ 
give the energy as a function of 5, 0, and a. For a 
given coupling strength g and angle Q, a{g, Q) is found 
by minimizing £. In general, this has to be done nu- 
merically. However, for small Fermi surface deformation 
a ~ 1 we can expand ([6|) in the small parameter t — v. 
The resulting integrals are straightforward and using 



J^^dv[U!{v)/v'~l]ir 



-0.9 we obtain 



1 + 1.065(1 - fcos^e) 
1 + 1.065(1 - |cos2e)' 



(7) 



Note that a < 1 as expected. 



Since £ 



it is now straightforward to evalu- 



ate the compressibility K{g,&). The resulting region 
of the collapse instability is shown in Fig. [2l There is 
a large region where the system is stable even though 
the interaction is strong (5 > 1) and partly attractive 
[8 < arccos(l/-\/3)]. This is most easily understood by 
ignoring the Fermi surface deformation {a — 1), which 
yields Qc arcsinl/v^ for 5 — > cx). In this limit, the 
interaction energy dominates and the gas collapses when 
the dipole-dipole interaction is attractive in more than 
half the x-y plane, i.e. for Q < arcsin l/-\/3. Fermi sur- 
face deformation effects are significant however; in Fig. [3] 
we show how a{g, Q) deviates substantially from 1 along 
the critical line Qdg) for collapse. 

The stability of the 2D system for a strong, attrac- 
tive interaction makes it a promising candidate to ob- 
serve superfluid pairing with unconventional symmetry. 
To investigate this, we solve the effective 2D BCS gap 
equation at T = 0, derived from the usual 3D BCS pair- 
ing Hamiltonian using V;(r) — 'i/;(p)$(z). This ansatz 
results in a 2D gap equation in terms of an effective 
interaction V2d{p) = J dz^'^{z)V{r), where ^r{z) = 
exp(— z2/2^^)/7r^/''Zy^ with 1^ — y/2jrnuj^ is the lowest 
harmonic oscillator wavefunction for the reduced mass 
m/2. Using the rescaled momentum k — (kx,ky) = 
{akx, ky/a) to describe pairing about the deformed Fermi 
surface, the 2D gap equation becomes 
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with 



k"^ /2m — /i, and /i the 



chemical potential. V2D(k, k') is the Fourier transform 

y2D(k,k') = y"dVsin(k.p)F2D(p)sin(k'.p). (9) 

of V2d{p), expressed in terms of the scaled momenta 
k,k'; i.e. V2d{K'^) = V'2n(k,k'). Note that only sin 
components contribute to the Fourier series since the or- 
der parameter is antisymmetric. The Pauli exclusion 
principle therefore cancels the short range r^^ diver- 
gence in the dipole-dipole interaction, making the Fourier 
transform fO]) finite even in the 2D limit. 

Using the expansions Aj^ ~ ^.^-^ A„(fc) cos[(2n — 
1)0] and l/2D(k,k') - E"„'=i K„' (fc, cos(2n - 
l)0cos(2n' — in ([5]), where is the angle between 
k and the x axis, we find numerically that Vnn ^ Vnn' 
for n' ^ n and all angles < < tt/2. This shows 
that to a very good approximation the gap is p-wave, 
Ak = Aocos0, and we can replace V2D(k, k') with 
Vii(fc, A;') COS0COS0'. Using this in ([H]), it reduces to 
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cPk' gs{k,k') cos^ (j)' 
(2^ 2E^, ■ 



(10) 



where Airgs/m = Vn. 

The dimensionless effective pairing interaction gs cor- 
responds to the dipole-dipole interaction averaged over 
the deformed Fermi surface, weighted by the p-wave sym- 
metry of the order parameter. After a straightforward 
but lengthy calculation, we find 



gs{k,k') = gM{k,k')Is[e,a{g,e)] 

9 

~ gM{k,k'){l - -cos^e). 



(11) 
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FIG. 3: The deformation parameter a and the coupling 
strength ga{kF,kF) for p-wave pairing along the critical line 
for collapse Qc{g)- For g < 0.725 where there is no collapse, 
we have evaluated a and gs at = 0. 



Here, Is and M are dimensionless functions given by 

/I I ...^^oJ,\/i 3 cos^ 0a^ cos^ \ 

d6{l + cos 20) (1 - ^, . .... J^, ) 
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and M{k,k') = 2171 [(1 + x)E{x) + (x - l)K{x)], for 
k < k' and k'lz <^ 1; for k' < k one simply swaps k and 
k'. K{x) and E{x) are the complete elliptic integrals of 
first and second kind respectively and x = k^ /k'^. We 
have M{kF,kF) = 1. For k'lz ^ 1, one has M{k,k') cx 
k{k'lz)~^ showing that the 2D effective interaction has a 
high energy cut-off for k'lz ^ 1- The last line in pT|) is 
exact for a = 1. 

We now analyze the T — gap equation pO)l as the 
dipole angle 8 is adjusted. Since M{k, k') > for all 
k, k', the gap Ak is zero unless /s(6, a) < 0. The result- 
ing superfluid region found by evaluating the sign of Ig 
is shown in Fig. [51 Numerically, we find that the Fermi 
surface deformation has little effect on the effective in- 
teraction for pairing as can be seen in Fig. [21 Ignoring 
Fermi surface deformation effects, we find from (jlip that 
the system is superfluid for 6 < 6s = arccos(2/3). The 
critical angle 8^ for superfluidity is therefore to a good 
approximation independent of the interaction strength g 
and is purely determined by geometry. 

Crucially, as 8 decreases from 7r/2, the system be- 
comes unstable toward pairing before the gas collapses, 
i.e. Qs{g) > ©0(5)) and there is a significant region in 
the phase diagram where the system is superfluid yet 
stable. Because the p-wave pair wavefunction is predom- 
inantly oriented along the x axis where the interaction 
is maximally attractive, the effective pairing interaction 
(ITT|) is stronger than the "bare" 2D interaction in ([5]) 
that determines the collapse instability. For this rea- 
son, even for strong interactions, there remains a window 
8s < 8 < 8c where the system is a stable superfluid. 

To examine the strength of the pairing interaction fur- 
ther, we make use of the fact that pairing occurs primar- 
ily at the Fermi surface, and in Fig. [31 we plot gs{kp, kp) 
from pT|) along the critical line Qc{g) for collapse. This 
gives the largest possible attractive pairing interaction, 
before the system collapses. We see that the effective 
pairing interaction increases monotonically with g and 
can become very large. Thus, one can produce a strongly 
paired gas without the system collapsing. This should be 
compared with the 3D trapped case where recent results 
indicate that the system is superfluid and stable only in 
a narrow region in phase space where the pairing is rel- 
atively weak [^, [13] ■ Solving the gap equation in the 
weak coupling regime yields Ao(8) = eF4e~^/2-\/Ae^/3= 
for 5s < (i.e. 8 < ©s) with A = k'^^^/kp an ultraviolet 
cut-off. This shows that the superfluid phase transition 



is inflnite order in the sense that 92 Aq] 



for all 



At flnite temperatures, long-range order is destroyed 
by phase fluctuations in 2D, and Ak(T > 0) = 0. 
The superfluid density remains finite, however, describ- 
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FIG. 4: Berezinskii-Kosterlitz-Thouless (Tbkt) and mean- 
field (T*) transition temperatures plotted as functions of O 
for g = 0.2. Inset: Temperature-dependence of the diagonal 
components of the superfluid density for = and g = 0.2. 



ing a Berezinskii-Kosterlitz-Thouless (BKT) superfluid 
of bound vortex-antivortex pairs [lH . The critical tem- 
perature Tbkt for this phase is the temperature at which 
the free energy of a single unbound vortex vanishes , 



T, 



BKT 



-Ps- 



(13) 



Here, ps = {ps,xx + Ps,yy)l'^ is the average of the diag- 
onal components of the superfluid mass density tensor 
Ps,ij- In estimating the energy E = \ J (Prps.ijViVj « 

In (■^) of a single vortex of radius a in a 2D box 
with sides of length L, we have assumed that the contri- 

due 



bution from the off-diagonal component ps, 
to the anisotropy of the vortex velocity field [13ll is small. 

lj| to allow 

for the effects of Fermi surface deformation, one finds 



Generalizing the usual expression for ps 



Ps,ij 



— mnS 



hkj- (14) 



As shown in the inset of Fig. [H rotational symmetry is 
broken by the external field and ps,xx > Ps,yy- The y 
component of the superfluid density tensor is suppressed 
since the quasiparticle spectrum is gapless in the direc- 
tion perpendicular to the field. 

Since ps depends on T, (fT3|) has to be solved self- 
consistently. We determine Ps{T) from (fT4|) using a T- 
dependent gap Ak(T) calculated by including the usual 
Fermi functions in the gap equation. The critical tem- 
perature determined this way is plotted in Fig. ([4]) for 
g = 0.2. We also plot the mean field transition temper- 
ature T* obtained from the gap equation. Phase fluc- 
tuations suppress Tbkt below T*. For weak coupling, 
Tbkt — T*. Since T* cx ^/ephajz and n2D ep, we have 



1 /2 

Tbkt cx . This suggests that one can cross the criti- 
cal temperature by adiabatically expanding the gas keep- 



^2D 



decreases (11 



ing T/ep constant so that T/Tbkt 
For stronger interactions, Tbkt quickly approaches its 
upper bound jmW^P ^ ^f/^Q oc n2D- 

Let us consider the experimental requirements for ob- 
serving the effects discussed in this letter. Recently, a gas 
of ^°K-*^Rb polar molecules was created with a density 
of n lO^^cm-^ at a temperature T/Tp > 2 [5]. These 
molecules have dipole moments of 0.57 Debye in their 
vibrational ground state. Writing g ~ 2(Pa~^NA/i00, 
where d = d/lDebye is the dipole moment in Debyes, 
a = a/lpm is the interparticle spacing measured in pm, 
and = m/niu is the mass of the dipoles in atomic 
mass units, the experimental parameters of Ref. 5] give 
g ~ 0.8. Also, the creation of 2D systems has already 
been achieved for bosons [l^. The observation of the 
effects discussed in this letter is therefore within experi- 
mental reach once further cooling has been achieved. 

In conclusion, we studied the quantum phases of a 
dipolar Fermi gas in 2D aligned by an external field. 
We demonstrated that by partially orienting the dipoles 
into the 2D plane, they can experience a strong p-wave 
pairing attraction without the system collapsing. This 
makes the system a promising candidate to study quan- 
tum phase transitions and pairing with unconventional 
symmetry. We also analyzed the BKT transition to the 
normal state at non-zero T. 
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